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Assuming the newly observed resonant structures X(1910), X(2200), and X(2350) as luoj, Locf), 

and (j>cj> molecular states respectively, we compute their mass values in the framework of QCD sum 

rules. The numerical results are 1.97 ± 0.17 GeV for uiu) state, 2.07 ± 0.21 GeV for uxj) state, and 

2.18±0.29 GeV for </>(/) state, which coincide with the experimental values of X(1910), X(2200), and 

X(2350), respectively. This supports the statement that X(1910), X(2200), and X(2350) could be 

^^ , uu), ucj), and (j>4> molecular candidates respectively. 

^H ■ 

Q ■ PACS numbers: 11.55.Hx, 12.38.Lg, 12.39.Mk 

>: 

O . I- INTRODUCTION 

:z;: 

f~^ " Very recently. Belle Collaboration reported observations of new resonant structures at M{lliijj) ^ 

CN ■ 1.91 GeV, M{lo4>) ~ 2.2 GeV, and M[(j)4>) ~ 2.35 GeV in 77 -> X ^ lou, w0, and (ficj) respectively 

,__i. [l|. For convenience, these resonances are called X(1910), X(2200), and X(2350) here. Since they all 

i-Q \ decay into two light vector mesons, it is natural to suppose them to be molecular bound states composed 

1^ of two light mesons. In theory, the molecular concept was put forward long ago ^2| and it was predicted 

^-1 that molecular states have a rich spectroscopy in Ref. [3[. The possible deuteron-like two-meson bound 

r~| ! states were studied in Ref. [J]. Although molecular states have not been confirmed in experiment, there 

'~~'" already have had some candidates for them. For instance, y(4260) could be a Yr/o" [5( or an uJXci state 



m 



X: 



Z+(4430) could be a D*Di state 0,0]; ^(3930) is proposed to be a D*D* [illH; y(4140) is interpreted 



> ; as a DlDl [13,0; ^(4350) could be a D^Df [li,Ei; y{4274)could be a DsD,o{2'in) 0. For 



more 



molecular candidates, one can also see some other Refs., e.g. 16|,ll7|. If molecular states can be completely 



o 

f — ■ confirmed by experiment, QCD will be further testified and then one will comprehend the QCD low-energy 

^*-^ . behaviors more deeply. Therefore, it is interesting to study whether the newly observed X states could be 

f^ \ molecular states. 

^vj ' In the real world, quarks are confined inside hadrons and the strong interaction dynamics of hadronic 

systems is governed by nonperturbative QCD effect completely. Many questions concerning dynamics of 
the quarks and gluons at large distances remain unanswered or, at most, understood only at a qualita- 



tive level. It is a great challenge to extract hadronic information quantitati vely from the rather simple 
Lagrangian of QCD. Fortunately, one can apply the QCD sum rule method [18| (for reviews see 19l422| 



and references therein), which is a nonperturbative formulation firmly based on QCD basic theory and has 
been successfully employed to some light four-quark states 23|-|29| . In Ref. [3C| , the authors have studied 



the tetraquark state qqqq by constructing and analyzing the sum rule composed of a diquark-antidiquark 
current with the quantum number J^ = 0+ and found masses of the tetraquark state qqqq appear in the 
region of 0.6 ~ 1 GeV, which are much lower than the mass of X(1910). Thereby, it may not likely to 
be a qqqq tetraquark state for X(1910). In Ref. [31j, the authors have studied the tetraquark udss of 
J^ = 0+ in the QCD sum rule and the mass of the tetraquark turns out to be around 1.5 GeV, which is 
much lower than the mass of X(2200). Thus, it may not likely to be a tetraquark udss for AT (2200). In 
Ref. 3^, the authors have predicted the mass of ssss tetraquark state of J^ = 0+ to be about 2.2 GeV in 
the relativistic quark model, which is slightly lower than the mass of X(2350). Just from the slight mass 
difference, one may not judge that A'(2350) is unlikely to be a ssss tetraquark state. However, one could 
at least see that the result does not exclude other possible interpretations such as molecular picture for 
^(2350). Therefore, we intend to obtain mass information of ujlo, uxj), and cfxj) bound states from QCD 



sum rules, and investigate whether ^^(1910), X(2200), and X(2350) could be new molecular candidates. 
The rest of the paper is organized as three parts. We discuss QCD sum rules for molecular states in Sec. 
im with the similar procedure as our previous works |33{ . The numerical analysis is made in Sec. IIIIl and 
masses of wcj, a;0, and (fx/) states are extracted out. The Sec. IIVI includes a brief summary and outlook. 



II. QCD SUM RULES FOR ojuj, u)4>, AND 00 MOLECULAR STATES 

The starting point of the QCD sum rule is to construct the interpolating current properly and then 
write down the correlator. In full QCD, the interpolating current for light vector meson can be found e.g. 
in Ref. [3J]. One can construct the molecular state current from meson-meson type of fields. Meanwhile, 
note that Belle Collaboration have indicated that there are substantial 0+ components in all three modes 
(77 — ?> X — ;> WW, ijjcf), and (t>(f))- Thus, following forms of currents with J^ — 0"*" are constructed for ww, 
ujcj), and 4>4> 

juu = {qci^qc){qc'it^qc'), (i) 

iu4, = {qci^qc){sc'itJ^Sc'), (2) 

34,4, = (Sc7^Sc)(Sc'7m^c'): (3) 

where q denotes light quarks u and d, with c and c' are color indices. One should note that meson molecules 
in the real world are long objects in which the quark pairs are far away from each other. The currents 
in this work and in most of the QCD sum rule works are local and the four field operators here act at 
the same space-time point. It is a limitation inherent in the QCD sum rule disposal of the hadrons since 
the bound states are not point particles in a rigorous manner. The two-point correlator is defined as 
n((?^) = i/d''a;e'''^(0|T[j(a;)j^(0)]|0). In phenomenology, the correlator can be expressed as 

. 2^ Af, 1 [^ , ImnP^^en/^) 
n(g ) = — T^ — o + - / ds 7^ + subtractions, (4) 

where Mh is the mass of the hadronic resonance, and \h gives the coupling of the current to the hadron 
{Q\j\H) = Xh- In the operator product expansion (OPE) side, the correlator can be written as 

,OPE 
s-q 

where the spectral density is p^ (s) — ^Imll^ (s), with the integration limit Smin ~ for ujuj state, 
Smin = (2?7is)^ for ujcf) statc, and Smin = (4?7is)^ for 4>(j) state. After equating the two sides, assuming 
quark-hadron duality, and making a Borel transform, the sum rule can be written as 

y2^^-MllM- ^ / ^^^OPEg-s/M^^^n^ond^ (g) 



"(?')-/ ds^— M+nCond(^2)^ (5) 



where M^ indicates the Borel parameter. To eliminate the hadronic coupling constant Ah, one reckons 
the ratio of derivative of the sum rule to itself, and then yields 

For the OPE calculations, we work at the leading order in as and consider condensates up to dimension 
ten, utilizing the light-quark propagator in the coordinate-space 



^^^ -^.~ ^^^.^ab9G,A^O^ +a'^^)~ -{qq} + — , 

- 3-^(39cr . Gq) + -^^r-^mg{gqa ■ Gq)^ ~ ^^q ^ (gg)(g^G^)- 



The s quark is dealt as a light one and the diagrams are considered up to the order mg. For some minor 
multi-gluon condensate contributions, one could omit them as the usual treatment. Concretely, spectral 
densities can be written as 

oTTCOnd ^ {gqcr-Gq){gqa-Gq) {qqfjg^G^} 
267r2 32.257r2 ' 

for ujuj state. 



p^'-Hs) = ^4^s\ P<-^\s)^^^s, P^^^^'is) ^ M_,^ ,(...-c.)(,)^M^:^ 



nisS, 



p{-){ff^G-^>(5) ^ _Mk^^^^ p{ss){gsa-Gs)^^^ ^ {s s) {gsa ■ G s) ^ ^(g«)(s«a-G?) (^) ^ {qq){gqa-Gq) 

X . VOtt'I 04™-z Q4/TT-Z 



^" ^ ^ — ^'^'^^ ^'"^ + 27^2 + 27,r2 + 32.267r2 + 32 • 2 V ' 

for ci;0 state, and 

oncond „ ^ ,5 ,n3 , {gsa-Gs){gsa-Gs) {ss)^{g^G^) 
mi - ms{ss) + 26^2 + 32.25^2 ' 

for (hd) state. 



III. NUMERICAL ANALYSIS AND DISCUSSIONS 

The sum rule ([7]) is numerically analyzed in this section. The input values are taken as nis — 



O.lOtaoa GeV [351, (qq) = -(0.23 ± 0.03)^ GeV^ {gqa ■ Gq) = ml (qq), (ss) = -(0.8 ±0.1) x (0.23 ± 
0.03)3 Qgv^ {gsa ■ Gs) = mg ^^^^^ ^2 ^ q g ^ g.l GeV^ and (52G2) = 0.88 GeV^ H- Complying 
with the criterion of sum rule analysis, the threshold -yio and Borel parameter Af 2 are varied to find the 
optimal stability window. In the standard QCD sum rule approach, one can analyse the convergence in the 
OPE side and the pole contribution dominance in the phenomenological side to determine the conventional 
Borel window: on one hand, the lower constraint for Af 2 is obtained by the consideration that the pertur- 
bative contribution should be larger than condensate contributions; on the other hand, the upper bound 
for Af 2 is obtained by the restriction that the pole contribution should be larger than the continuum state 
contributions. Meanwhile, the threshold ^^so is not arbitrary but characterizes the beginning of continuum 
states. For many hadrons, the first excitation of studied state defines the size of i/io, and the difference 
between ^/s^ and the mass Mh of studied state is around 0.5 GeV. Concretely, the value of ^/sq is fixed 
by these steps: 1) taking a value of ^/sq; 2) fixing the corresponding Borel parameters Ar2 according to 
two rules (OPE convergence and pole dominance); 3) extracting the mass from the sum rule in the work 
window fixed in the first and second steps; 4) checking that whether the ^/sq chosen in the first step is 
acceptable using the empirical relation that the difference between ^/sq and Mh is around 0.5 GeV; 5) if 



the ^/sq chosen in the first step is not acceptable, return to the first step, vary ^/sq and go on. Taking 
UJLU as an example, we choose ^/so = 2.4 GeV and finally arrive at Mjj = 1.97 GeV. One could check 
that y/so = 2.4 GeV is acceptable with the empirical relation that the difference between ^/sq and Mh is 
around 0.5 GeV. Thus, we choose the central value of y^ = 2.4 GeV for ujlu state. 

Here, we would make some particular discussions on the choice of Borel windows. It has been shown 
in detail in some Refs, such as |30] and [36|, that it is not possible to find an conventional Borel window 
for some light scalar tetraquarks. The problem is that the four-quark condensate is very large, making 
the standard OPE convergence to happen only at very large values of AP. In fact, it has appeared the 
same problem in this work. Taking the ujuj as an example, the comparison between pole and continuum 
contributions from sum rule (J6|) for ujuj state for ^/sq = 2.4 GeV is shown in the left panel of FIG. 1, and 
its OPE convergence by comparing the perturbative with other condensate contributions is shown in the 
right panel. Even if we choose some uncritical convergence criteria, e.g. the perturbative contribution 
should be at least bigger than each condensate contribution, there is no standard OPE convergence up 
to M^ > 1.8 GeV^. The consequence is that it is unable to find the conventional Borel window where 
both the OPE converges well (i.e. the perturbative contribution bigger than each condensate contribution) 
and the pole dominates over the continuum (the latter one happens at M^ < 1.3 GeV^). Under such a 
circumstance, one could try several possible ways to solve the problem. I) One could release the criterion 
of pole dominating over continuum and take some high values of Borel parameter M^. Thus, OPE series 
can converge well. Graphically from the Borel curve, one can see that there is a very stable plateau. Some 
authors have virtually adopted this way to deal with the above problem. However, there occurs some 
other problem. Although there is very good OPE convergence and a flat plateau for the Borel curve, 
contributions from continuum states are dominating. As one knows, the phenomcnological side of the sum 

rule can be expressed as n(g^) = j^'i"_ 2 + ^ /j°° ds *"" ^_ ^ due to the "single-pole+continuum states" 
hypothesis. From the criterion of pole dominating over continuum, one can obtain the maximal value of 
the Borel parameter M^ satisfying the "single-pole+continuum states" model. Exceeding this value of Af ^ , 
the single-pole dominance condition will be spoiled. Thereby, the Borel parameter A'P must not be chosen 
too large to warrant pole dominance. II) One could push the threshold parameter ^/sq to a very large 
value, and the maximum value of M^ will be enhanced with the increasing of ^/sq. Thus, one may find the 
Borel window satisfying both the perturbative bigger than condensate contributions and the pole bigger 
than continuum contributions. However, the threshold parameter y^io is not arbitrary but characterizes 
the beginning of the continuum states. With too large values of ^/sq, contributions from high resonance 
states and continuum states may be included in the pole contribution. Hence, the QCD sum rule may not 
work normally. HI) One could warrant the pole dominance firstly and try releasing the strict convergence 
criterion of perturbative contribution larger than each condensate contribution in some case. In the present 
work, we have dealt with the problem in this way. It is worth to note that the treatment is not arbitrary 
but there is some definite condition. For example, we consider the ratio of perturbative contribution to 
the "total OPE contribution" (the sum of perturbative and other condensate contributions calculated) 
but not the ratio of perturbative contribution to each condensate contribution. Not too bad, there are 
two main condensate contributions with different signs (four-quark condensate and two-quark multiply 
mixed condensate) and they could cancel with each other to some extent, which brings that the ratio of 
perturbative contribution to the "total OPE contribution" is bigger than 60% at M^ > 0.8 GeV for ujuj for 
y/sQ = 2.4 GeV. In addition, we calculate and find that the ratio of perturbative contribution to the "total 
OPE contribution" does not change much including some high dimension condensate contributions. In this 
sense, we could say that the OPE converges in the region while satisfying pole dominance. Although it may 
not be good OPE convergence in comparison with the conventional case, one could find a comparatively 
reasonable work window. Note that the treatment should not be arbitrarily transplanted to any case. One 
could take the Ref. [36] as an example. From its FIG. 4, one can see that for /o all condensate contributions 
in the region 0.4 < M^ < 0.7 GeV are larger than the perturbative contribution. That means one could 



not even find a region that the perturbative dominates in the "total OPE" allowed by the upper bound. In 
a word, to deal with the problem on choosing the conventional Borel window in QCD sum rules, which have 
similarly appeared in some other multiquark states, we warrant the pole dominance preferably and release 
the strict convergence criterion to a weak one that perturbative dominates in "total OPE contribution" , 
so that the convergence of OPE is still under control while satisfying pole dominance. Although it may 
not be so good OPE convergence as the conventional case, one could find a comparatively reasonable work 
window and extract the hadronic information of studied states reliably. Thus, we choose the minimum 
value of M^ to be 0.8 GeV^ and the maximum AP to be 1.3 GeV^ for ww state for ^5^ = 2.4 GeV. 
Similarly, the maximum value of M^ is taken as 1.2 GeV^ for y/so = 2.3 GeV; for y/so = 2.5 GeV, the 
maximum M^ is taken as 1.4 GeV . The dependence on M^ for the mass of ujoj state from sum rule ^ 
is shown in FIG. 2, and we arrive at 1.97 ± 0.07 GeV for cjuj state. Considering the uncertainty rooting 
in the variation of quark masses and condensates, we gain 1.97±0.07±0.10 GeV (the first error reflects 
the uncertainty due to variation of y^ and Af^, and the second error resulted from the variation of QCD 
parameters) or 1.97 ± 0.17 GeV for ujuj state. 

The comparison between pole and continuum contributions from sum rule ([6|) for uxjj state for y/Jo — 

2.6 GeV is shown in the left panel of FIG. 3, and its OPE convergence by comparing the perturbative 
with other condensate contributions is shown in the right panel. There has the same problem for ujcj) as 
the above case for uiui, and we treat it similarly. For ujcj) state, the ratio of perturbative to the "total 
OPE contribution" at M^ = 1.1 GeV^ for ,/s^ = 2.6 GeV is around 67% and increases with the M^. 
Furthermore, the relative pole contribution is approximate to 54% at AP — 1.5 GeV^ and descends with 
the M^. Thus, the range of A'P is taken as AP = 1.1 ^^ 1.5 GeV^ for ^/sq = 2.6 GeV. Similarly, the proper 
range of M^ is obtained as 1.1 ^ 1.4 GeV^ for y/s^ = 2.5 GeV, and the range of M^ is 1.1 -- 1.6 GeV^ 
for ,/so — 2.7 GeV. The mass of ujcfi state as a function of M^ from sum rule ([7]) is shown in FIG. 4, 
and we obtain 2.07 ± 0.13 GeV for cucj). Varying input values of quark masses and condensates, we attain 

2.07 ± 0.13 ± 0.08 GeV (the first error refiects the uncertainty due to variation of ^/sq and AP, and the 
second error resulted from the variation of QCD parameters) or 2.07 ± 0.21 GeV for cucf) state. 

For (j)(j} state, the comparison between pole and continuum contributions from sum rule © for ^/sq ~ 
2.7 GeV is shown as an example in the left panel of FIG. 5, and its OPE convergence by comparing the 
perturbative with other condensate contributions is shown in the right panel. A bit difference for the 
case of 00 is that the perturbative contribution can be bigger than the second most important condensate 
(^sf' at AP > 1.1 GeV^. Meanwhile, the pole contribution can dominate in the total contribution 
while AP < 1.6 GeV . Thus, it is possible to find a region where both the OPE can converge well 
(the perturbative contribution bigger than each condensate contribution) and the pole dominates over the 
continuum. Thus, the range of Af ^ for 00 state is taken as AP = 1.1 ^ 1.6 GeV^ for ^/Jq = 2.7 GeV. Via 
the similar analyzing process, the proper range of M^ is obtained as 1.1 ^ 1.7 GeV^ for y/so = 2.8 GeV, and 
the range of M^ is 1.1 ~ 1.8 GeV^ for y^ = 2.9 GeV. In the chosen region, the corresponding Borel curve 
to determine the mass of 00 state is shown in FIG. 6, and we extract the mass value 2. 18 ±0.20 GeV for 00 
state. Subsequently, we vary the quark masses as well as condensates and arrive at 2.18 ± 0.20 ± 0.09 GeV 
(the first error refiects the uncertainty due to variation of ^/sq and A4'^ , and the second error resulted from 
the variation of QCD parameters) or 2.18 ± 0.29 GeV in a concise form. 

IV. SUMMARY AND OUTLOOK 

In 77 — > A — > LUio, aj0, and 00, Belle Collaboration observed three new resonant structures at A/I{LUio) ^ 
1.91 GeV, A'I{lu(J)) ^ 2.2 GeV, and A/(00) ^ 2.35 GeV. Assuming these newly observed resonances as 
molecular states, we have employed the QCD sum rule method to calculate their masses, taking into account 
contributions of operators up to dimension ten in the OPE. Our final numerical results are 1.97±0.17 GeV 
for WW state, 2.07 ±0.21 GeV for w0 state, and 2.18 ±0.29 GeV for 00 state, which are in agreement with 





FIG. 1: In the left panel, the solid line shows the relative pole contribution (the pole contribution divided by the 
total, pole plus continuum contribution) and the dashed line shows the relative continuum contribution from sum 
rule © for ^/so = 2.4 GeV for uju state. The OPE convergence is shown by comparing the perturbative with other 
condensate contributions from sum rule ((6| for ^/so — 2.4 GeV for lolj state in the right panel. 
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FIG. 2: The mass oi ujuj state as a function of M from sum rule ((7| is shown. The continuum thresholds are taken 
as ^ = 2.3 ~ 2.5 GeV. For ^ = 2.3 GeV, the range of M^ is 0.8 ~ 1.2 GeV^; for 
M^ is 0.8 ~ 1.3 GeV^ for ^ = 2.5 GeV, the range of M^ is 0.8 ~ 1.4 GeV^ 
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FIG. 3: In the left panel, the solid line shows the relative pole contribution (the pole contribution divided by the 
total, pole plus continuum contribution) and the dashed line shows the relative continuum contribution from sum 
rule ((6l for y^ = 2.6 GeV for lo(I> state. The OPE convergence is shown by comparing the perturbative with other 
condensate contributions from sum rule ((6]) for y^ = 2.6 GeV for aj(j!) state in the right panel. 
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FIG. 6: The mass of 00 state as a function of M^ from sum rule ^ is shown. The continuum thresholds are taken 
as yi^ = 2.7 ~ 2.9 GeV. For ^i^ = 2.7 GeV, the range of M^ is 1.1 ~ 1.6 GeV^; for ^ = 2.8 GeV, the range of 
M^ is 1.1 ~ 1.7 GeV^; for ^) = 2.9 GeV, the range of M^ is 1.1 ~ 1.8 GeV^ 
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the experimental values of X(1910), X(2200), and X(2350) respectively. This supports the statement that 
X(1910), X(2200), and X(2350) could be cjw, w0, and (fxp molecular states respectively. However, one 
should note that there are still some differences between our central values and experimental values. At 
present, we have merely considered ujuj, uxj), and (fxf) molecular states with J^ = 0+. Belle Collaboration 
indicated that while there are substantial spin-0 components in all three modes (namely 77 — )• X — > ujuj, 
uj(f), and (fxj)), there are also spin-2 components near threshold. The differences between our central values 
and experimental data are probably caused by that we have not considered the spin-2 components for 
UJUJ, uj<j), and (jxj) state here, which implies that the theoretical predictions might be improved by including 
J = 2 components for the future. In addition, one needs to take into account other dynamical analysis to 
identify the nature structures of these X States for further work. 
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